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Abstract 

Remote state preparation is generation of a desired state by a remote observer. In spite of causality, it 
is well known, according to the Reeh-Schlieder theorem, that it is possible for relativistic quantum field 
theories, and a “physical” process achieving this task, involving superoscillatory functions, has recently 
been introduced. In this work we deal with non-relativistic fields, and show that remote state preparation 
is also possible for them, hence obtaining a Reeh-Schlieder-like result for general fields. Interestingly, in 
the nonrelativistic case, the process may rely on completely different resources than the ones used in the 
relativistic case. 
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I. INTRODUCTION 


Remote State Preparation (RSP) QH3) is a generation of a desired state by a remote observer, 
who merely acts locally on his own distant system. Such a process is considered successful, if the 
remote observer is able to make sure that for particular choices of a measurement and its outcome 
the remote system is in the desired state. Although the success probability might, in general, be 
small for a single run of the process, it is required that for events with a successful measurement 
result, the remote state will be arbitrarily close to the desired one (i.e., with fidelity arbitrarily close 
to 1). 

A causal theory is a theory with a maximal velocity c, due to which information cannot be 
transported instantaneously across a distance R, but rather at some finite time T > R/c. Causality is 
an important property of relativistic systems, and in particular of relativistic quantum field theories 
(which are usually referred to only as ‘quantum field theories’ - QFTs). As the standard model of 
high energy physics is a quantum field theory, the importance of causality in our universe and the 
physical theories describing its behavior could not be overestimated; However, when one wishes 
to relate relativistic theories, and in particular QFTs, with RSP, he could, naively, argue that the 
concept of RSP is impossible in a causal theory. 

Surprisingly or not, according to an important theorem established long ago by Reeh and 
Schlieder [4}-;6] in the context of algebraic quantum field theory (AQFT), RSP is possible in rela¬ 
tivistic quantum fields. According to the theorem, the action of a field operator, or a polynomial of 
it, on the vacuum is dense in the Hilbert space - i.e., allowing to generate states which are arbitrar¬ 
ily close to any desired state in the Hilbert space, regardless of its location in configuration space. 
Recently 0, a physical process resulting with RSP of any desired state was presented, suggesting 
a “physical” way to derive the theorem (rather than AQFT): a prescription for the remote prepara¬ 
tion of a desired state was given, rather than a general proof of its existence. This prescription is 
based on superoscillatory functions, suggesting a deep relation between the task of RSP in QFTs 
and the mathematical phenomenon of superoscillations [® 0- As the fields are relativistic, it is 
clear that the resource which allows for such a task is entanglement [|T0HTT| . For other quantum 
fields, which, in general, do not necessarily posses an entangled vacuum state, it remains unclear 
whether such a task could be accomplished. 

One might ask then what is the case for nonrelativistic fields. On one hand, these may be 
noncausal, and hence one shall not be surprised of noncausal effects such as remote state prepara- 
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tion. On the other hand, how can one be sure that this still allows to remotely prepare any state? 
In this work we show that this task is possible for a very large class of quantum field theories 
which satisfy some reasonable physical assumptions. It turns out that RSP is possible, even for 
fields with a non-entangled vacuum and a bounded group velocity. Thus, a Reeh-Schlieder-like 
result is obtainable for general quantum fields. However, for fields with a non-entangled vacuum 
state, the mechanism which produces the remote states relies on an infinite front velocity as a re¬ 
source, rather than entanglement. This has interesting physical consequences. On one hand, the 
entanglement-based mechanism, allows generating field states even with a presence of a “wall” 
between the operating region and the target region, as nothing actually propagates in the process. 
This is not possible in case the mechanism is based on infinite front velocity. On the other hand, 
the latter mechanism is less sensitive to noises caused by previous attempts to perform the process. 

This paper is organized as follows: First, we shall formulate the physical problem we address 
and describe the process of state preparation in section II. In section III we argue when superoscil¬ 
lations are required and discuss their implications on the success probability. Finally, in section 
IV, we make the connection to the Reeh-Schlieder theorem and argue that the processes described 
hereby is a Reeh-Schlieder-like result in the sense that it allows generation of arbitrary field states, 
even though it may utilise completely different physical resources. 

II. STATEMENT OF THE PROBLEM 

Consider a quantum field theory in d+l dimensions, whose action is invariant under translations 
and rotations (reflections for d = 1), but not necessarily under boosts. The dispersion relation of 
such a theory, due to the rotational invariance, satisfies 

to = to(k), (1) 

with k = |k|. Due to the translational invariance, the eigenfunctions take the form / k (x, /) = 
h (to) e ,(k ' x ~ a,n and thus the field operator may be expanded as 

0(x,f)= f (to k ) (a k e i<k ^ n + a k e- Hk ^), (2) 

using annihilation and creation operators, satisfying the canonical commutation relation [ a k , a q | = 
(2 n) d 6 {d) (k - q). We also assume that h (co k ) is a real function - otherwise, its phase may always 
be absorbed in the creation and annihilation operators. 
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The field operator ([2]) is the solution of the Heisenberg equation of motion derived using the 
field hamiltonian H {) . For our purposes, we couple the field to a detector - a two level system, 
located in x = x 0 , described by the Hamiltonian H d = Q.cr z /2 (fi = 1). The interaction between the 
field and the detector is given by 


H int = A (cr+e (t) + cr_e* (t)) <p (x 0 ), 


(3) 


and is only switched on for a finite period of time - e (t) =£ 0 only for t e [-/ 0 ,0]. We assume that 
in t < -to the field is in its vacuum state, and the detector - in its ground state, ||). Assuming that A 
is very small and that e{t) is 0(1), one may calculate the state in t = 0, after the interaction, using 
first order interaction picture. There 


H int = 2ie (t) cr + e lQ,t (p (x 0 , t) + h.c.. 


and thus the detector-field state after the interaction is given by 


(4) 


l<M> = |0. 


,1 )-a f 

J-tc 


dte it) e iQt (p (x„, t) |0, T). 


(5) 


Post-selecting the spin in the excited state IT), one obtains the (normalized) field state 

d d x 


|®> 


= «'*f 


:e(co k + Q)h(to k )e ,k x |k), 


(2 n) 1 


( 6 ) 


where e(aj k + £2) = JT° dte(t) e i(bJk+cl)t , M = f h 2 (co k ) \e(co k + f2)| 2 and |k) = a{ |0). 


-to J J (Inf 

Every single excitation of the field may be written as 


m = n: 


f 


d d xF (x) (p (x) |0). 


(7) 


where AA = f d d xd d y-^h 2 (oj k ) F* (y) e ,k (y X) F (x) [12.]. 

J (Z7r) 

We wish to choose a window function e(t), such that the state IT 7 ) will be created as a result of 
the interaction. Thus, we demand that |( V F|®)| = 1. This is satisfied if we choose 


^des (&A- + f2) 


/■ 


d d xe~ ik x F (x), 


( 8 ) 


where ‘des’ stands for ‘desired’; i.e., the Fourier transform of the desired window function is equal 
to the one of F (x). 

The spherical symmetry restricts us to the creation of states which are symmetric around the 
detector. Thus, F(x) = F (r), where r = |x - x 0 |; WLOG, we assume next that x 0 = 0. In 
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this case, the integration reduces to 6des(<4t + £2) = f drr d ~' Fir) f dO suT* 2 Q e ~ lkrcose J d d ~ 2 Cl. 
After performing the angular integrals f d d ~ 2 Q = ( r u |/+ \)p } (f° r derivation see llT3l ) and 
f ddsin d ~ 2 Q e ~ ,krcose = (kr I2) {2 ~ d)l2 T ((J - 1) /2) /«- 2)/2 (kr) (where J v (x ) is a Bessel func¬ 
tion), one obtains 


?des (.OJk + fi) 


f 


drkF ( r) 


2 nr\ 


d/2 


Jd^g (kr). 


(9) 


This equation can be satisfied if and only if a> k = a> ( k ) is bijective. As r increases, the oscillations 
of the Bessel function (in k space) become faster, hence e des (a> k + Q) obtains larger and larger 
Fourier components. At some point, say r > r 0 , the function e des (a> k + Q) begins to have significant 
Fourier components which oscillate (in frequency space) faster than t 0 (below we show that r 0 is 
related to the minimal group velocity of the theory). Since e (/) is nonvanishing only within [-/ 0 ,0], 
the standard frequency-time relations of Fourier transforms suggest that the above relation could 
not be satisfied for r > r 0 . 

It was shown in 0 that, for relativistic fields, this problem can be circumvented using a super- 
oscillatory e(co k + Q). In the following section we introduce this solution briefly. Then, we show 
that while the same solution holds for non-relativistic fields, some consequences of it are different. 
This allows us to answer the question “when are superoscillations required?” in a general manner. 


III. SUPEROSCILLATIONS 

Superoscillatory functions are functions that oscillate faster than their fastest Fourier compo¬ 
nent BH |9j. This is due to a destructive interference, and thus they are always accompanied by 
exponentially larger amplitudes somewhere outside the so-called superoscillatory region. In our 
implementation, we would have to place the exponentially larger amplitudes in a non-physical do¬ 
main of e. Another difficulty regarding superoscillations is that these functions can superoscillate 
in an arbitrarily large (but not infinite) domain. Therefore, there must also be a physical non- 
superoscillatory domain. In our implementation, we shall choose a superoscillatory function that 
will not be exponentially amplified in this non-superoscillatory domain. Then, we would have to 
find a way to eliminate the contribution of the function in this domain. 

We shall now proceed by finding functions which meet these demands. Consider the following 
function lf71 [T4l IT5ll : 

~e m W) = — r dae^^e* cos( “"' A) , (10) 

26 \2n Jo 
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where D, 8 and A are some constants, and to' = to + Q-too for some ojq. Since t = t 0 (cos a - 1) /2, 
the function e |/!| (t) has support only in [-/ 0 , ()]. Nevertheless, we shall now prove that e [/l] (to ) 
oscillates arbitrarily fast in to' e [0, m c ] for to c <sc (s 2 t 0 cosh[A]) . Performing the integration 
explicitly we obtain 


e [h] (to') = 


V2d 


e~^ t0 Jn 


1 

d 2 


1 +8 2 to't 0 cosh [A] + -8 4 to' 2 tl 


(ID 


For to' > 0. Using the asymptotic form of the Bessel function [fl6l for 8 <s: 1, and then taking 
6 2 «: (6u c t 0 eosh [A]) -1 we obtain 


~e lh] (to') = De 


.11 /r 

' C ° S * S 2 + 2^ C ° Sh ^ ~ 4 


( 12 ) 


in this domain. Redefining e [/,] (to') to be the summation of two such functions, one having 8 2 = 
2 7tm + jt/4, and the other D —> +/D and d -2 = 2yr/n - ^/4, where m » 1, we get 


gW (oj) = £)e^ iw t ’o( ±c °sh[A]-i) 


(13) 


This function oscillates in to space at “frequency” t = (t {) (± cosh [A] - 1). By increasing A we 
can set these oscillations to be arbitrarily fast. The superoscillatory domain is finite, therefore the 
condition described in Eq. 0 cannot be exactly satisfied. However, one can get arbitrarily close to 
satisfying this condition by increasing the superoscillatory domain. This is achieved by decreasing 
8 . 

Superoscillations come at the price of an exponential growth outside the superoscillatory do¬ 
main. In our case the growth occurs at to' < 0. For fields whose energy is bounded from below, 
we choose to 0 such that to > to 0 , therefore to' < 0 corresponds to to + Q < to 0 , which is in the 
non-physical domain. Beyond the superoscillatory domain, the function gradually obtains regular 
(slower) oscillations, and in the limit to' » to c , it behaves like o/ 1/2 sin (to'to). In order to elimi¬ 
nate this contribution, one can convolute (t) with some function h (t) which is differentiable n 
times (n » 1) and has a small temporal support. 

Finally, we use a combination of such superoscillatory functions, each with a different t', in 
order to generate the window function 

e(to') = J' dt'e lh] (to';t')e des (t'). (14) 

In the limits T —» oo and 8 —> 0 we get e(to') —> e des (to') in the segment to' e [0, to c ]. (This is 
while the actual window function, e (/), and the desired window function, e des (t), are very different: 
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e(t ) has temporal support only in [-/ 0 , 0 ], while 6 des (t) might have an arbitrarily large temporal 
support.) Therefore, we can generate remote spherical symmetrical one-particle field states around 
the spin, up to an arbitrarily small infidelity. The generalization to arbitrary field states is achieved 
using the same method introduced for relativistic field states 0. Here we only note that the 
generalization to one-particle states which are not spherical symmetrical involves an array of 
spins (located at different positions in an arbitrarily small region) rather than a single spin, and the 
generalization to many-particle states involves a set of such arrays. 

It is shown in 0 that the exponentially small amplitude of e (co') in the superoscillatory domain 
results in a success probability of the form 

u c T 2 

P~e~^r. (15) 


The finiteness of the superoscillatory domain is responsible for an infidelity 


n 



(16) 


Since F (k) is normalizable, in the limit > oo one obtains 77 —» 0. Inverting the latter functional 
relation to = u> c (t]) = 1 /g(rf), and expressing T as a functional of the desired state |'F), we get 
the relation 


P ~ g sW' 0 


(17) 


When F (k) decays with a power law, g ( 77 ) behaves according to a power law as well, and when 
F (k) decays exponentially g(rj) ~ 1 /In (1 ///). For a given dispersion relation, 7’ [| V F)| can be 
calculated by plugging e des (oJ k + Q) = dte‘ {0H:+ - l)1 (0 in the l.h.s of Eq. (9) and expressing 

the r.h.s as a Fourier transform of a temporal function [00. 

When are superoscillations required? Mathematically, the answer can be deduced from Eq. 
0 . However, in order to understand the meaning of that, let us first consider a very simple, 
non-superoscillatory, window function, e(t) = 6 (t + to) - a short and impulsive interaction, for a 
quantum field theory in 1 + 1 dimensions. Assuming the detector is in the origin, the generated 
state of the field, after the detector’s post-selection, will be 


X oo 

dkh (u) k ) (a[ k e~ i0Jkt ° + | 0 >. (18) 


Both terms represent wave packets propagating out of the detector: one is left-moving and the 
other is right-moving. Assuming that co k 6 R Vk, and that the group velocity is less than the phase 
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velocity, i.e., dco/dk < oj/k, each wave with k, co reaches at t 0 distance 



( 19 ) 


(the minimal r 0 (co) is the r 0 mentioned above in the simple case of 1 + 1 dimensions and a specific 


desired state). Thus, only waves with frequencies satisfying v g (co) t 0 > L propagate fast enough to 
arrive to ±L without superoscillations. In relativistic theories, as well as other theories in which 
the group velocity is bounded from above, one can always find L such that all frequencies would 
require superoscillations, while in theories in which the group velocity is unbounded from above, 
for every arbitrary L , some frequencies will not require superoscillations. 

The above case deals only with waves which are outgoing from the detector. However, one 
may also consider the case of ingoing wave packets - for example, if one wishes to generate the 
field state 


ITO = (cf>(L) + cp(-L)) |0) 



( 20 ) 


which contains both ingoing and outgoing wavepackets. The ingoing waves require v g (co) < 0 for 
k > 0, which is impossible, and thus superoscillations are also required for this case. 

For general desired states, L is roughly the separation between the operating region and the 
farthest place in the target region. Note that the same description holds for higher dimensions; in 
these cases (when using a single detector) the resulting state is spherically symmetrical. While 
every single point on the sphere generates waves which propagate in all directions, the spherically 
symmetrical state as a whole generates waves which propagate only in the +r directions due to the 
Huygens principle. Thus, superoscillations are required for the generation of wavepackets which 
propagate faster than the (slowest) group velocity, or inwards, into the interaction region. 

IV. RELATION TO THE REEH-SCHLIEDER THEOREM 

It is now the right time to recall a long-standing, possibly surprising result of AQFT - the Reeh- 
Schlieder theorem 0-161. According to this theorem, the set of Hilbert space vectors \\J/) e 'H 
generated from the vacuum (or any other bounded state with a bounded energy BH) of a relativistic 
QFT by operating with polynomials of the field operators in any open region is dense in 'H. In 
other words, by applying certain local operators to the vacuum state in a certain region O i, one is 
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able to generate, with nonzero success probability, a state of the field localized at some remote re¬ 
gion^) {O k } k > 2 , arbitrarily close to some desired state. The {O k } k > 2 regions may remain, throughout 
the process, outside the light-cone of 0\, and thus this outcome must be the result of pre-existing 
vacuum correlations (the theorem entails a violation of Bell inequalities as well lfT8l [1911. Con¬ 
sider, as an example, a relativistic scalar field <p with mass m. There, for |x - x'| 2 » m~ 2 , one 
obtains the equal-time correlation function 

(0| (p (x, t) (p (x', t) |0) ~ e~ m ^ x ~ x '\ (21) 

- the correlations do not vanish even between spacelike separated regions [|20ll . 

In this paper, we have shown that RSP is theoretically possible for general quantum fields, 
including non relativistic ones. Thus, we have described a Reeh-Schlieder-like process for general 
fields. Note that while the Reeh-Schlieder theorem does not involve time dependence (as one 
would expect in a relativistic context), our process does. In relativistic theories, adding time 
dependence is meaningless, as for any two spacelike separated points in spacetime, x /J = (/, x) and 
x' IJ = (/', x'), there is a reference frame (connected by a Lorentz transformation), in which t = t' 
and x?tx'. Thus, if we define AjC = x? - x' M , and r 2 - -Ax M Ax^ = - (t - t') 2 + (x - x') 2 > 0, due 
to Lorentz invariance, in the limit r 2 » m 2 , one obtains (0| (p (x, t) (p (x', /') |0) ~ e~ mr . In the other 
extreme case of fields which do not possess any correlations at all - i.e., fields for which 

<0| cp (x, t) <p (x', t) |0) oc 6 {d) (x - x') , (22) 

the time dependence is crucial, because correlations are generated in time. For example, the 
Schrodinger field satisfies 

(0| (p (x, t) (p (x', t) |0> = d W) (x - x') , (23) 

and thus, without the time dependence, the overlap between the generated state and the desired 
state is zero and RSP is not possible. Therefore, one can deduce that RSP is possible only when 

(0| (p (x, t) (p (x\ t') |0) ^ 0, (24) 

for every x, x’ and t ± t'. It means that for every quantum field (in which RSP is possible) either 
there exist correlations in t = t' and/or some components of the quantum field propagate infinitely 
fast, i.e., the front velocity is infinite. This is yet another formulation of a result discovered in 
mm and widely discussed in [23] [24]. In the first scenario, our mechanism uses vacuum 
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correlations as a ‘resource’ for RSP, while in the second scenario the ’resource’ is the infinite front 
velocity. Remarkably, in both cases the (gedanken) prescription for RSP is the same. 

The different ‘resources’ have different physical implications. Consider for example a case 
where one puts a wall between the operating region and the target region. A good model for such 
a wall, in the case of a scalar field, for example, is a time dependent potential, 

// wa „ = limT f i d xW (x) 6 (t + t 0 ) (p 2 (x), (25) 

r-»o° J 

where 6 (t + to) is the Heaviside step function and W (x) = 1 where the wall is placed and 0 
elsewhere. This potential adds an “infinite mass” to the field in particular space points where 
W(x) = 1, starting from t = -t 0 , and thus it acts as a wall. 

In causal theories with the equation of motion 0(p (x IJ ) = 0, the Green’s function defined by 

OG - /) = -id {4) (x M - /), (26) 

may be chosen to be a causal, retarded Green’s function, G (x^ - y M ) = D R (x M - y M ), which van¬ 
ishes when y ,J is outside the past light-cone of x 1 ', and thus the field state during the preparation 
process can be expressed as a superposition of field states at time t = -t 0 - all in the past light-cone 
of the generating region. Therefore, if the wall is placed outside the past light-cone of both the 
operating region and the target region throughout the process, its effect, which is propagating at 
the speed of light at most, will not travel far enough to destroy the vacuum correlations between 
the two regions and RSP would be possible. 

On the other hand, when the process is due to infinite front velocity, the wave front will en¬ 
counter the wall and RSP will not be possible. Mathematically, this is manifested in the fact that 
in this case of a non-causal field, there are no light-cone and “causal” Green’s function. Therefore, 
the field state cp (x, -t 0 + e) |0) (where e > 0 is arbitrarily small) will involve field states at / = -t 0 
from regions where the wall is placed. 

While the previous example was in favor of the vacuum correlations of relativistic fields, one 
can also come up with different examples which favor fields with an infinite front velocity. Con¬ 
sider, for example, a case in which one fails to remotely generate a field state and then tries again. 
The unsuccessful attempt contaminates the vacuum state of the field. Due to high order terms in 
Eq. ([6]), in relativistic fields, the field state will no longer be of bounded energy (i.e. it will be a 
superposition of energy states with ever increasing energies) ll251 . This will destroy the delicate 
vacuum correlations and so RSP will not be possible. However, for some non-relativistic fields, at 
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FIG. 1: A wall which is placed at -to between the operating region (denoted as O i) and the target region 
(denoted as Oi) will disrupt the process for fields with an infinite front velocity. However, for relativistic 
fields, if the wall is placed outside the past light-cone of both 0\ and Oi throughout the interaction, the 
process will not be affected by it. 

least, one could repeat the process on and on. This can be easily demonstrated using fields with a 
bounded frequency. 

In both scenarios, as correlations exist for every x, x' pair, superoscillations should not be in¬ 
terpreted as the generators of overlap between the initial state and a desired, remote one. Super¬ 
oscillations should be rather understood as cancellations of the correlation of the initial state with 
any other one, except for with the desired state. 

Since remote preparation of field states can be used to generate entanglement between distant 
observers, this result implies that every quantum field 

1. which is invariant under rotations and translation, 

2. whose spectrum is bounded from below, 

3. whose dispersion relation is bijective, and 

4. whose front velocity is bounded from above 

- possesses a vacuum state which violates Bell’s inequalities. 
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V. DISCUSSION 


In this paper, we have discussed remote state preparation for non-relativistic fields, and argued 
that it was possible. This was explicitly shown using a physical process. The scheme that we have 
considered in this paper is similar to the scheme discussed for relativistic fields in 0. It involves 
coupling the field to a spin (or an array of spins) and then postselecting the spin to a certain state 
which corresponds to the creation of one (or more) quanta of the field. 

In relativistic fields, due to vacuum correlations, particles are never completely localized. I.e., 
particles created in any arbitrary point in space have a nonvanishing overlap with particles created 
anywhere else in space. In nonrelativistic fields, although particles may be localized when created, 
after an infinitesimal time they have a nonvanishing overlap with particles created anywhere else 
in space as well. I.e., localized particles have an infinite front velocity. As both vacuum correla¬ 
tions and infinite front velocity induce particles generated in the “operating region” to have some 
(usually very small) overlap with particles generated in the “target region”, both are mechanisms 
upon which the process of field RSP may rely. 

In order to cancel the overlap between the initial state and states in any other location, except 
for the target region, a specially tailored window function is used. When the target region is 
too far for the wave components with the slowest group velocity to reach during the interaction 
duration, the window function must be superoscillatory. The farther the target region is, the faster 
the superoscillations have to be. This holds for any theory, either relativistic or not. The success 
probability and fidelity of the process depend only on the frequency of the superoscillations, and 
thus, in particular, they do not depend on the specific mechanism or the existence of a light-cone, 
but rather on the dispersion relations of the field. The specific mechanism and the existence of a 
light-cone become important when the process is slightly altered. For example, when placing a 
wall between the operating region and the target region before the process begins. 
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